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We present analytical perturbative, and numerical solutions of the 
Einstein equation which describe a black hole with a nontrivial dila- 
^ ■ ton field and a purely topological gauge potential. The gauge potential 

has zero field strength and hence no stress-energy, but it does couple 
to virtual string worldsheets which wrap around the Euclidean horizon 
two-sphere, and generate an effective interaction in the spacetime la- 
grangian. We use the lagrangian with a nonstandard potential for a 
scalar field that reproduces the effect of the worldsheet instantons. As 
has been previously pointed out the topological charge Q of the gauge 
q-i field can be detected by an Anharonov-Bohm type experiment using 

quantum strings; the classical scalar hair of the solutions here is a classi- 
cal detection of Q. ADM mass, dilaton charge and Hawking temperature 
are calculated and compared with the known cases when appropriate. 
We discuss why these solutions do not violate the no-hair theorems and 
show that, for sufficiently small interaction coupling, the solutions are 

stable under linear time dependent perturbations. 
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1. Introduction 

The no-hair theorems in the physics of black holes essentially state that the only physical 
parameters describing a black hole are its mass, angular momentum, and conserved charges 
associated with long range fields. Originally no-hair theorems arose from the uniqueness 
theorems of Israel, Carter and Wald [ttl for the Kerr solution. If the cosmic censorship 
conjecture is correct, and hence the gravitational collapse of a star gives rise to a black hole, 
the no-hair theorems imply that the black hole must be described by the Kerr solution. 
An alternative method of proving a no-hair theorem, developed by Bekenstein || is to 
show that the black hole geometry cannot support nontrivial classical field configurations 
that are of interest. We review this in section 5. In particular, a no-hair theorem must 
be proved (or disproved) for each matter theory of interest. For a scalar field coupled 
to gravity, there are simple proofs when either the potential energy V((p) or <fidV/d(j) is 
positive definite polynomial. So scalar hair is ruled out for the standard self interacting 
potential V(<f>) = m 2 <p 2 /2 + \<f> /4 if m 2 , A > 0. The important case of spontaneous 
symmetry breaking, when m 2 < and <j> is coupled to a gauge field, requires a more 
careful arguement, and is first proved in ||. 

In recent years, however, black holes with various kinds of hair have been discovered. 
The "colored" black hole supporting Yang-Mills gauge field was first numerically discovered 
0, and existence was later proven analytically || . The colored black hole possesses zero 
Yang-Mills charge and thus parametrized only by its ADM mass. Further, it was shown 
|| that a black hole can exist inside a magnetic monople of spontaneously broken gauge 
theories. Under much more generality it was shown in that black holes can exist 
inside various classical field configurations of which magnetic monopole is only a special 
case. The charged dilaton black holes || || carry electric and magnetic charges, and also 
have nontrivial scalar fields. Now, in all these examples of black holes with scalar hair, 
the black hole also carries electric and/or magnetic charges. In this paper, we present 
numerical solutions, and analytic perturbative solutions, for black holes with nontrivial 
dilaton field, but no electric or magnetic charge. Instead, the black hole has a topological 
"axion" charge, that is, there is a nontrivial 2-form gauge potential B a b, with vanishing 
field strength, that wraps around the horizon two-sphere. Classically B a b is not observable; 
its existence can be inferred by a distant observer, by the nontrivial scalar field on the black 
hole. 

Thes e are asymptotically flat versions of the black holes constructed by Giddings 
et.al. [Kj and are the solutions to bosonic string theory, which correspond to black holes 
with a culaton and a purely topolgical axion gauge potential. However, these exact solutions 
are "all throat," i.e. the radius of two sphere is constant. The authors write down a 
spacetime lagrangian which contains an effective potential giving the effect of Euclidean 
string world sheets that couple to the gauge potential and wrap around the horizon two- 
sphere. It is this interaction term which allows one to get around the usual no hair 
construction. Hence the classically observable scalar hair is supported by the classically 
unobservable world sheet instantons. 

In classical electrodynamics the gauge invariant field strength is the physically ob- 
servable quantity. However in quantum mechanics the potential can affect a particle even 
in those regions where the field strength is identically zero, as illustrated in the famous 
Aharonov-Bohm effect. A similar phenomenon comes up in string theory. A string natu- 
rally couples to a 2-form potential, B a b, since its trajectory is described by a world sheet. 
H a bc(= V ' a Bbc + cyclic terms) is the gauge invariant quantity and determines the classical 
equation of motion. The static spherically symmetric solution of an axion field interact- 
ing with gravity is the Schwartzschild metric supplemented by an axion potential with 
a nonzero component given by [[ll| Bq^ = a'Qsm.6/2 where Q is a constant. It follows 
that H = and no classical measurement outside the black hole could determine the axion 
charge of a black hole. However as pointed out in |nj if two strings originating at the same 



spacetime point are made to interfere at another spacetime point with the volume enclosed 
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by two world sheets enclosing the black hole, the quantum phase difference, between two 
strings at the end of their trajectory will be proportional to Q. Thus Q becomes a new 
observable quantity of the system despite the fact that strings world sheet do not touch 
the region of nonzero H. This observable could be considered a new kind of quantum hair 
for black holes. 

The Anharonov-Bohm type experiment described above for strings would be difficult 
to perform. If the axion field were coupled to another auxiliary field, such as the dilaton, 
one might infer axion charge by measuring the auxiliary field. Interestingly one can give 
such an interpretation to an effective spacetime action descri bing the interaction between 
gravity, the axion and the dilaton derived by Giddings et.al.[|10|. In their work an exact 



classical solution of bosonic string theory is found as a conformaTly invariant sigma model, 
containing both the symmetric metric coupling, and an antisymmetric coupling Bg^ = 

a' Qsm O/2. A fixed point of that theory occurs when Q = n. These exact solutions found 
in |IIJ can be interpreted as a black hole geometry without an asymptotically flat region. 
As such their solution cannot describe a black hole found in our spacetime. However, 
on large scales the system can be studied using a spacetime action with the world-sheet 



instanton effects included by an effective interaction [1C]. In the rest of the paper we 



investigate the existence of asymptotically flat solutions using this spacetime action. 

2. Equations of Motion and Boundary Conditions 

We will write the static spherically symmetric Einstein metric as 

dv 
ds 2 = -f(r)dt 2 + ^— + R 2 (r)dn 2 (2.1) 



The four-dimensional low-energy effective action obtained in [|10| for Q = n is 



S = ( 4^j [-11 + 2(V<pf + V ((/>)] d 4 x 



where <p is the dilaton field and V reproduces the effect of the the world sheet instantons 
which couple to the topological gauge potential Bg^, 

Vtf) = 2 _9_ e ^-2R\r)l a ' (22) 



Note that we are using the Einstein metric, whereas in |T(| the string metric is used, 
dimensionally reduced to the t-r plane. 



In (|2.2|) , C is a positive constant which comes from the determinant of fluctuations 
about the instantons. Hence C naturally arises as a small parameter, and this suggests 
using a perturbative approach to look for solutions which are almost Schwarzchild black 
holes, with a small but nonconstant dilaton. We first find these solutions, and then turn 
to the exact case. 

After choosing the gauge R = r, we obtain three independent equations of motion, 
with r = xrh 

Ax 2 <P" = (1 + x(f)')Ce 2 ^- x2r ^ - (1 + A)x<f>' (2.3) 

A'x = l- A- Ax 2 (/)' 2 - Ce 2 ^- x2r ^ (2.4) 

(In/)' = 2x(P' 2 + (lnAy (2.5) 



where primes denote differentiation with respect to x, and 2/a' has been scaied away. For 
asymptotically flat spacetime, we require 

A, f -»• 1 - 2M/r as r -> oo. (2.6) 

Requiring an event horizon at r = r^ (a; = 1) implies that 

A = h = (2.7) 

The subscript h will stand for quantities at x — 1. For a solution we require that </> is difler- 

entiable at Tk. The energy density, p, is given by 4ivp = —4tvT} = (Ac/)' 2 + Ce 2 ^~ r /r 2 )/2. 
Finiteness of p at r^ puts no further condition on <p. The total energy in the dilaton field 

outside the horizon is finite, Air J pr 2 a/ '/ /Adr < oo, if <p' vanishes asymptotically at least 

as fast as 1/r with 5 > 1.5. By solving ( |2.3|) asymptotically, we find that 5 = 2. Thus 

<P ->• 0oo + 0(l/r) as r -> oo (2.8) 

together with equation (|2.6|) constitute sufficient conditions for aysmptotic flatness and 
finite total energy. To get a unique solution we must impose on a particular value of </> at 
infinity. We can require that 

0oo=O (2.9) 

though we will see in section 4 that if C is fixed there is a range of <ph and hence a r ang e 
of 0OO that leads to the asymptotically flat solutions. Alternatively we can impose ( |2.9| ), 
without loss of generality, for all solutions and let C vary in which case there will be a 
range of values of C that give ( |2.9| ) . That there is one-to-one correspondence between the 
two parameters, (p^ and (7, can be easily seen from the fact that setting C = 1 amounts 

to redefining <p according to the scheme Ce 2 ^ — > e 2< ^ . Each of the two parameters that 
one allows to vary has its own advantage. If one keeps (poo fixed and lets C vary, one can 
imagine turing on the perturbation potential V off the Schwartzschild black hole. But to 
discuss the existence of solution numerically it is more convenient to fix C and let (p^ (or 
equivalently <ph) vary 

3. Perturbative Solutions 

We observe first that when C = the equations of motion for the metric functions are 
solved by the Schwartzschild spacetime. The equation for (p is that of a massless scalar 
field minimally coupled to gravity and is solved by the unique solution </> = constant. For 
small C we look for a perturbative solution to the equations of motion. Non-existence of 
a perturbative solution would indicate that no exact solution exists while existence means 
only that the exact solution may exist. 

Let the metric be that of Schwartzschild with f = A = 1 — 1/x and C and to be 

small quantities of the same order. The equation of motion for (p, 2V 2 = V, after scaling 
r = xrh, becomes 

[(x 2 - x)4>']' = Ce 2(j} - x2 ^ 
which can be immediately integrated to obtain <p to first order in the small quantity 

<p=c r^idt 

Joot 2 -t 

where G(z) = ^e~ Th ^ d(. This solution satisfies all boundary conditions: -£■ ~ 1/x 2 , 

<fi ~ as r — > oo, and 6 is bounded and differentiable at r%. It is also consistent with our 
assumption that 4> an d C are of the same order in the small quantity. 

One can also compute the first order corr ectio n to the metric. If one lets 
A = 1 — (1 + C8)/x then one can solve for 5 using ( J2.4J) . One finds 

5 = G(x) 

Since lim^^op G(x) = positive constant, this solution indicates that the presence of the 
dilaton field increases the ADM mass from the Schwartzschild v alue , a fact supported by 
the numerical solution found in the next section. From equation (|2.5| ) it follows that / = A 
to first order in C. 



4. Numerical Results 



For all cases considered in this section, we take rh — C = 1 unless otherwise specified. 
For all other values of r^ we have tested, we find that numerical solutions also exist, but 
we do not have a proof that this is true in general. At the end of this section we discuss 
the parameter space of fixing (p^ and varying C. 

Because equatio ns (|2.3| ) and ( |2.4| ) do not involv e / w e sol ve th ese two equations for A 
and 4> and then use fl2.q) to solve for /. Equations ( |2.3fl and (|2.4[ ) are solved numerically 
with the boundary conditions ( |2.6| ), and (|2.7|) using the shooting method |12]]. To start 



shooting, initial values must be chosen for A, 
gives a constraint between <ph and cj)' h : 



and <j>' at rh- But at x — 1, equation (|27 



<t>'h = V(e r 



-2<j} h 



-1) 



(4.1) 



Note that the solution in section 2 satisfies this constraint. For a given value of r^ we thus 
have only one shooting parameter, <fih- Equation ( |4.1| ) is singular when 2<ph = r\. Indeed, 
a solution exists only when 2<ph < r\, because the equation ( |2.4j ) gives 



A' h = l-e 2< ^- 



(4.2) 



and since Ah = it is necessary that A' h > which is possible only if 2(ph > r\. 

One sees that <\> monotically increases with x as follows. For 2<j)h > r^, 4>' h > 0. Thus 
if there exists a local extremum of 6 for some x = x e > 1 it must be a local maximum. 
However according to equation (|2.3|) 

A( Xe )x 2 e <f)"(x e ) = e mx«)-* 2 y h > 



i.e. that local extremum can only be a local minimum. Therefore there can be no local 
extremum of (p fo r any finite value of x. 
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Figure 1. A vs. x 
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(r h ,C) = (1,1) 
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Solid line is the Schwartzschild solution, 1 — \JX. For the dashed lines, from 

left to right, (p h = -0.4, 0.1, 0.2, 0.235. 



We find that for each r^ there is a range of 4>h £ ( — °°? <^™ ax ] that will give 
^oo =constant. For example for r^ = 1, <ph = —0.21678 gives (p^ = and </>™ ax = 0.235596. 
The relation between <^>™ ax and r^ is summarized in figure 7. For all <ph < ^ iax <p monot- 
ically increases to a constant value and A approaches the Schwartzschild result, 1 — 1/x, 
asymptotically. Figure 1 also shows that A approaches Schwartzschild solution everywhere 
as <ph — ► — oo. This is expected because 4>h — > — oo is the same as C — > + . 

Once A and are known / can be integrated using equation ( |2.5| ) . Note that (|2.5| ) is 
linear in /. This is compatible with the fact that the time coordinate can be scaled by a 
constant. Because of the event horizon, (|2.5|) does not determine f' h . Initially we take f' h = 
A' h which gives a solution such that linv^oo / = /oo. Further scaling of time coordinate 
is necessary to obtain linv^oo / = 1 i.e. for the final solution we pick f' h 
necessary for computing Hawking temperature later. 
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Figure 2. <fi vs. x 



For the numerical analysis it is convenient to think of the solutions as being 
parametrized by (rh,<f>h)- However, we may also consider the solutions as being 
parametrized by two global charges: the ADM mass, M, and the dilaton charge, D. We 
now turn to calculating these two charges. 

The ADM mass can be read off of the large r behavior of A, 1 — 2M/r. Figure 3 shows 
that for a fixed </>oo, (M/Mc=o) increases with C as noted with the perturbative solution. 
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Figures. (M/M c =o) vs. C (r h , <^o) = (1,0) 



A plot of (x 2 <p') vs. x in figure 4 shows that <p' ~ constant/x 2 as x — > oo. Thi s 
asymptotic behavior agrees with that of perturbative solution and satisfies equation ( |2.8| ) . 
With this behavior of (j)' , 

is finite. D should not be thought of as another free parameter. Figure 5 shows one-to-one 
correspondence between D and (p^. Also there is one-to-one correspondence between two 
parameter spaces (r h,(f> oo) and (M,D). 
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Figure 4. (x 2 (f)') vs. x (r h , C, (j)^) = (1, 1, 0) 
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Figure 5. D 



Another quantity of interest is Hawking temperature, Th = k/2-k. k is surface gravity 
, j md is given by \im r ^ rh (Va) where a = ^a c a c , a c = (£ 6 V b £ c )/(-£ d £ d ), £ c = (-§1)°, 
and V = y/-€ c £ c - We find that 



A-kTh = lim f'/y/fjA = Jf[# h = A'J^JZ 

continuously decreases with (p^, and hence with (ph- An alternative expression, obtained 
by using equation (|4.2|), 4ttTh = (1 — e 2 ^ >h ~ rh )/y / f 00 explains the decreasing of Th with 



(f>h- It also supports our finding that no solution exists for large value of <f>h because if 
4>h > 4>c = r ^/2, Tjj will be negative. Figure 7 shows the relation between <^>™ ax and r^ 
and it shows that as r^ — > oo, ^>? iax — ► ^^/2 = Hawking temperture bound. 

For a fixed r^ the effect of dilaton charge is to lower the Hawking temperature as 
evident by figure 6. In this respect the dilaton charge has a similar effect on the black hole 
as the electric or magnetic charge. Recall that the Hawking temperature of a Reissner- 
Nordstrom black hole is given by 



4ttT, 



Q 



i 

em 



H 



where r_|_ is the outer horizon radius, and Q e 
hole. 



the electric or magnetic charge of the black 
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Figure 6. AixTh vs. D r\ l = \ 
Throughout this section we have taken C = 1 and found that the solution exists when 
i>h G (— oo, </>™ ax ] (or equivalently (j)^ G (— oo, 0™ ax ]-) As noted at the end of section 2 if 

i> is the solution with C = 1 then </> will be the solution for a different positive C where 

i = (j) — In yC. If we demand that (j)^ = always, then for a fixed r^ the solution exists 



if and only if C G [0, C max ] where C max = & 
d>h — 4>h — <$>oo and C = e 2 ^°° will ensure 6^ - 



It can be easily verified that picking 



0. 



5. Consistency with No-Hair Theorems 

To check whether black hole solutions we found are consistent with various no-hair the- 
orems we review the proofs such theorems. Consider a theory of scalar field interact- 
ing with gravity described by the action, S = f[— 1Z + (V</>) 2 /2 + V(<j), x^yf—gd^x. If 
we write the line element in the form of (|2.1|) with R = r the equation of motion is 

\/A/f (r 2 (f)' \J f A) = r 2 dV/d(j) where primes denote d/dr. Now for asymptotic flatness we 
can consistently require linv^oo <p = 0. Multiplying both sides by <p an d integrating by 



parts gives 



¥r 2 y/JA 



oc 



r r*VTZ 4> n + %+/ 'A 



dr. The left side is zero on account 

dv 



of the boundary conditions imposed at r^ and oo. If ^-0 is positive definite for r > r^, 
then the integrand is positive definite and we see that for such forms of V((f>) the only 
solution for </> is the trivial one: = 0. However, for the interaction V we use (p^fr is not 
necessarily positive definite. 



No-hair theorems using scaling arguments [13] impose slightly different condition, 
namely, that V itself be positive defi nite . It appears at first sight that V we use satisfies 
this condition. However the proof in [13] makes use of the fact that V does not depend on 
coordinates explicitly which our V clearly does not satisfy. 



6. Stability Analysis 

To study the stability of the solution, we will consider spherically symmetric perturbations 
around the static solution. We willl show that the solution is stable for sufficiently small 
C. Let 

f{r,t) = f {r) + 8f(r,t) 

A(r,t) = A (r) + 6A(r,t) 

6(r,t) = 6 (r) + 86(r,t) 

Q(r,t) =TT + 6Q(r,t). 

/o, Aq, and 6q are the static solution found in the previous sections and 8 denotes smallness 
compared to static solution. Here we are making a perturbation around the topological 
charge Q as well, which in the static case was taken to be n. The new action is 



S 



- n + 2 (V</>) 2 + e" 4 ^ - ^ e ^-2rVa' CQsQ 



-gd^x 



where H a b c has nonzero components 



H r 



/Shi# 
a o r 6Q 



and 



H 



te<j> 



a ——d t 8Q 



The linearized gravitational field equations obtained from the action do not lead to 
an unstable mode because a spherically symmetric gravitational field has no dynamical 



degrees of freedom. Explicit demonstration of this fact is given in \TA . 

After dropping the subscript of the background solution, thelmearized equation of 
motion for 86 is, 




2 2 



./ 



x 2 



where ' = d/dx (x = r/rt) and a' has been scaled away. We have also scaled the time 
coordinate by a constant and taken the time dependence of 86 to be e ut . Then, after 
transforming to the tortoise coordinate -^ = ^- r , we obtain the equation for ip(y) = xl ' 

as 

d 2 tjj 



where E 



-UJ 



and 



fA- 



dy' 



+ (E- V^tP = 



V* 



ix ax x^ 



24>-x 2 r 2 h 



This is one dimensional nonrelativistic schrodinger equation for a particle with mass = | 
moving under the influence of the potential V^. Instability of the solution corresponds to 
the existence of the bound states, E < 0, solutions of the schrodinger equation. But V& is 
positive definite for most of the parameter space of the static solution. See figure 7. Thus 
in this region of parameter space, where C and/or r^ are small, a linear time dependent 
perturbation in the dilaton field leads to no instability. 1 
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Figure 7. 4>h vs. r^ 



According to the remarks at the end of section 5 this figure can be looked upon 
as O vs. T/j with (poo — 0- Dashed line is the curve 7Y/2, Hawking temper- 
ature bound on (ph with Cy = 1. Solid line is (p^ vs. T\i\ it determines the 
boundary of solution space i.e. no static solution exists for any point above 
the solid line. Vd, is positive definite below the small-dashed line. 

We next investigate the equation of motion for 5Q in the parameter space in which V<f, 
is positive definite. Assuming the same time dependence as before the linearized equation 
of motion for 5Q can also be put into the form of the schrodinger equation as before with 
the potential energy function 



V Q (y)=Af[B'< 



B'(x)]-^(fA)-3fx 2 rf l Ce^- 



where B(x) = 2cp'(x) + 1/x. Here we have set a' = 2 for convenience. To study Vq it 
is convenient to think 2 of the representation in which (r/j, C) is varied and (p^ = 0. One 
should imagine figure 7 as C vs. rn with cph —*■ —oo being C — > 0. Now the last term in Vq 
will be small when C and/or r^ are small whereas other terms do not depend explicitly 
on C or r/,. Thus the effect of the last term in Vq will show up for the points in upper 
right of figure 7. Figure 8 shows a typical Vq when C and/or r^ is small; it has a well - 



We believe that even in the region of parameter space which makes V<j> negative no bound 
state with E < is possible. The reason is that when V^ does become negative it is in the form 

of a shallow well surrounded by a high barrier and a small barrier. 

Because Ce 6 ^ appears in Vq there is no simple tranformation between the parameter space 
where C is fixed and ^ varied and that where ^ is fixed and C varied. However, when one 
considers which part of the static solution space is stable under perturbation, it is consistent to 

use the simple transformation rule given at the end of section 5. 



9 



which we will refer to as the 
consider varying (rh,<j>h)- 



'primary" well — followed by a barrier. From now on we will 




Figure 8. Vq vs. x (r h , C, <f> h ) = (0.077161, 1, -1) 
The effect of r^ and 4>h when they are larger (i.e. when they are in the upper right of 
figure 7) is to introduce a "secondary" well whose location depends on r/j. In general this 
secondary well gets deeper as r^ and (f>h gets bigger and moves to the left as r^ gets bigger. 
See figures llabcd. If there is any bound state of Vq the particle will be localized in the 
primary well or the secondary well. Since big r/j and <ph give a deep secondary well, we 
start looking for bound states for such cases. 

To find the bound states of Vq it is numerically simpler to use the original x coordinate. 
As x — > oo the solution for 5Q is e~ ux and the linearized equation of q = 5Qe LJX is 



q" - (a + 2u)q' + (u 2 + uoa - b)q = 



where 



d/dx, 



2B(x) 



l^'(fA) 



and 









The wave function in the schrodinger equation is related to q via ip = (q/x)e~ 2 ^~ UJX . 

For arbitrary value of to and arbitrary potential Vq, no regular solution of ip will be 
possible. We call ifj a regular solution if it goes to zero sufficiently fast as x — > 1 and 
x — > oo. We find that for a deep enough secondary well of Vq and a discrete set of positive 
values of to regular solutions exist. These regular solutions are the bound states we are 
looking for. 




-0.06 



100 1000 

(r h ,C,<f> h ) = (0.077161,1,-0.875) 

Solid line is the wave function of the ground state and horizontal line the energy level. 
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Figure 9. ip and Vq vs. x 
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As shown in figure 9 the wave function is localized around the secondary well, and this is 
evidence that no bound state will be possible when the secondary well is shallow enough. 
Therefore for a fixed r^ there may be bound states when <ph is close to 0™ ax but as 
(/>™ ax — ► — oo bound states will disappear. For a fixed r^ we can graph u> of the ground 
state against 4>h and by extrapolating find the zero of that graph, 0°. See figure 10. For 



h< 



no bound state is possible and hence the solution is stable. 



7. Summary 

We have presented the evidence for the existence of scalar hair of a black hole in the 
presence of string instantons which couple to a topological gauge potential. The space time 
action that we use includes an effective interaction due to the string instantons wrapping 
around the Euclidean black hole. Measurement of the dilaton charge of the black hole can 
be looked upon as an indirect detection of the topological axion charge of the black hole. 
For example, if a "dilatonic ball" held at a large distance from the black hole experiences 
an attractive force towards the black hole in excess of the gravitational force the dilatonic 
charge of the black hole will be ascertained, and hence the axionic charge. 
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Figure 10. to vs. 



(r h ,C) = (0.077161,1) 



b° h « -0.892 



Because there is one-to-one correspondence between G and <ph (of the static 
solution) this figure can be roughly looked upon as 00 vs. O with (poo — 0. 
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Figure 11a 



(0.0005265,1,-1.33594) 



Figure lib 




(0.003679,1,-1.194) 



Figure lie. Vq vs. x (r h , C, 4> h ) = (1, 1, 0.23) 
i>/j's of figures a,b,c lie on the small dashed line of figure 7. 




Figure lid. Vq vs. x (r h , C, <j> h ) = (1, 1, -1.194) 
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